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The archetype of geometrically frustrated compounds SrCrgpGai2-gpOi9 is a kagome bi-layer 
of Heisenberg Cr'^"'' ions (5 = 3/2) with antiferromagnetic interactions. We present an extensive 
gallium NMR study over a broad Cr-concentration range (0.72 < p < 0.95). This allows us to 
probe locally the susceptibility of the kagome bi-layer and separate the intrinsic properties due to 
geometric frustration from those related to site dilution. Compared to the partial study on one 
sample, p = 0.90, presented in Phys. Rev. Lett. 85, 3496 (2000), we perform here a refined study 
of the evolution of all the magnetic properties with dilution, with a large emphasis on the lowest 
diluted p = 0.95 sample synthesized for this study. 

Our major findings are: 1) The intrinsic kagome bi-layer susceptibility reaches a maximum at 
a temperature of « 40 — 50 K, which we show here to be robust to a dilution as high as « 20%. 
This maximum is the signature of the development of short range antiferromagnetic correlations in 
the kagome bi-layer ; 2) At low-T, a highly dynamical state induces a strong wipe-out of the NMR 
intensity, regardless of dilution; 3) The low-T upturn of the macroscopic susceptibility is associated 
to paramagnetic defects which stem from the dilution of the kagome bi-layer. The low-T analysis of 
the p = 0.95 NMR lineshape coupled with a more accurate determination of the nuclear hamiltonian 
at high-T, allows us to discuss in detail the nature of the defect. Our analysis suggests that the 
defect can be associated with a staggered spin-response to the vacancies of the kagome bi-layer. 
This, altogether with the maximum in the kagome bi-layer susceptibility, is very similar to what 
is observed in most low-dimensional antiferromagnetic correlated systems, even those with a short 
spin-spin correlation length; 4) The spin glass-like freezing observed at Tj ~ 2 — 4 K is not driven 
by the dilution-induced defects. 



PACS numbers: 75.30.Cr, 75.50.Lk, 76.60.-k 

I. INTRODUCTION 

Under certain circumstances, it is impossible for mag- 
netic systems to minimize simultaneously all the interac- 
tions between the spins. The system is then frustrated. 
This is known to occur in the spin glass (SG) compounds 
where the disorder on the magnetic network induces a 
competition between the interactions. In a vast variety 
of systems, the frustration can arise,from the geometry of 
the lattice itself, without disorder ,EJ as in the case of the 
triangular-based antiferromagnetic (AF) networks. The 
kagome {d = 2) and the pyrochlore {d — 3) networks 
with AF interactions are a particular class of geomet- 
rically frustrated networks where the triangles (for the 
kagome lattice) and the tetrahedras (for the pyrochlore 
lattice) share corners instead of sides as for the familiar 
triangular network. 



Within a classical theory, the ground state is built on 
triangles and/or tetrahedras with a zero total magnetic 
moment. Remarkably, the corner-sharing geometry giyegi 
rise to a macroscopic degeneracy of the ground statefl'EI 
By macroscopic, we mean that the energy of the kagome 
or of the pyrochlore ground state is invariant under a 
rotation of a finite number of spins. Through these ro- 
tations, the system can explore all the spin configura- 
tions making up the ground energy level. The energy 
spectrum is therefore characterized by the existence of 
zero energy excitations, the so-called soft modes. □ In par- 
ticular, the soft modes associated to a small number of 
spins are extremely efficient in destroying a long range 
magnetic order. A non magnetic ground state is pre- 
dicted at T — > K, with a spin-spin correlation func- 
tion (SiSj) ^ exp{—rij/S,), where the correlation length 
^ does not exceed twice the lattice parameter. El The huge 
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reservoir of soft modes leads also to a low energy shift 
of the excitation spectrum, so that an wiusually highly 
dynamical state is predicted at low-TB This peculiar 
ground state is commonly called a cooperative param- 
agnet state or a "spin liquid state" . 

Quite remarkably, a different description of quantum 
nature for S* = 1/2 geometrically frustrated AF networks 
yields similar results. As first elaborated by Anderson,Ll 
and validated through numerical calculations, the ground 
state can be constructed on spin singlet pairs, and not 
on non magnetic triangles (or tetrahedras) .□ In par- 
ticular, numerical studies point at the existence of a 
"gap" in the magaetic excitation spectrum of the kagome 
S = 1/2 petworkjj and maybe in the pyrochlore S* = 1/2 
network. t3 The term gap is actually inappropriate as a 
continuum of singlet states is imbedded between the sin- 
glet ground state and the first excited triplet state. These 
singlet excitations are in some regards the analog of the 
classical soft modes. A spin liquid ground state is again 
favored, with unusual thewnodynamic properties such as 
a high entropy at low-T.llil 

The discovery of the pkagome-based insulator 
SrCr9pGai2_9pOi9 in 1988,E3 brought considerable 
attention on geometric frustration. Since then, an 
intense mapping of the low-T physics of the kagomej-and 
of the pyrochlore compounds has been carried outJlj'til 
Original properties haiie been uncovered, encompassiiag 
anomalous SG states|l3 and ice-like ground states.t£l 
However, nowadays, only a few compounds are good 
candidates for a spin liquid ground state. Following Vil- 
lain's early work,llZl the discrepancy between experience 
and theory is to be found in the perturbations to the 
ideal AF Heisenberg hamiltonian with nearest neighbpi: 
spin-spin interactions. The dilution of the netwock,EIE£l 
the interactions 0|ther than nearesi- neighbor,t2l the 
dipolar interactionj^j the anisotropyJEj etc. are pertur- 
bations which differentiate the geometrically frustrated 
compounds and make them deviate from the ideal 
spin liquid behavior. Each deviation can potentially 
induce a long range order, usually quite complex. In all 
experimental studies, one must therefore discriminate 
between the properties related to geometric frustration 
and those related to these perturbations, that we may 
label by the general term of disorder. 

However, even in the presence of these limiting param- 
eters, spin liquid-like compounds do exist .E3 More than in 
any other compound, spin liquid-like properties are ob- 
served in SrCr9pGai2-9pOi9 [SCGO(p),0 < p < 1]. The 
geometric frustration in SCGO arises from a kagome bi- 
layer of Heisenberg Cr'^^ ions {S=3/2), a quasi— d = 2 
network of two kagome layers connected by a triangu- 
lar lattice linking layer (Fig. |^). The disorder in SCGO 
stems from the dilution of the Cr-network by non mag- 
netic Ga'^^ ions. To date, all SCGO crystals reported in 
literature are non stoichiometric {p < 1). 

Although an anomaly is observed at low-T {Tg w 2 — 4 
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FIG. 1: The magnetic lattice of SCGO is a stacking of kagome 
bi-layers of Cr"^"'" ions (« 78% of the Cr sites), separated by 
Cr-Cr isolated spin pairs (~ 22% of the Cr sites). 



experimental data available on SCGO point at the ex- 
istence of a spin liquid-like ground state. The neutron 
diffraction pattern at T < Tg is characterized by a broad 
peak, from which is extracted a spin-spin correlation 
length of twice the Cr-Cr distance (^ « 2dcr-Cr)!^ A 
more refined neutron experiment showed that the diffrac- 
tion pattern is in agreement with the existence of sub- 
groups of spins (singlets, triangles or tetrahedras) of zero 
magnetic moment EH Neutrons also revealed that only 
a fraction of the Cr'^+ moment is frozen and /zSR ex- 
periments revealed [tie existence of a strongly fluctuat- 
ing ground state,B3'E3 a picture further supported by re- 
cent specific heat measurements that indicate that puiy 
^ 50% of the total entropy is removed below 100 K.E3 

From the macroscopic susceptibility {Xmacro)^ o^i^ can 
extract that a strong AF interaction couples the neigh- 
boring spins with a characteristic Curie- Weiss temper- 
ature Qmacro ~ 500 — 600 K. The Curie- Weiss behav- 
ior of Xmacro SxtCnds tO tcmperaturCS T ^ Qmacro, 

which is remgirized as the most typical signature of 
frustration.E3o Puzzlingly, at very low-T, Xmacro P^'O" 
gressively deviates from Curie- Weiss behavior and ex- 
hibits a behavior closer to a simple Curie law. This prop- 
erty of SCGO is actually encountered in the majority of 
the geometrically frustrated magnets. This is in deep 
contrast with the non magnetic ground state predicted 
for a spin liquid. It was conjectured that this Curie up- 
turn might originate from the dilution, i.e. may not be aa 
intrinsic property of the kagome bi-layer susceptibility.Ell 
If this is indeed the case, what is then the susceptibil- 
ity of the kagome bi-layer? And what is the underlying 
mechanism which triggers the Curie upturn? 

These are the main topics addressed in this work 
through the NMR of gallium nuclei of SCGO. The ^^Ga 
and ''^Ga nuclei (/ = 3/2) are local probes coupled to 
the Cr'^^ ions. The gallium nuclei labeled Ga(4/) are at 
the heart of frustrated physics, as they are exclusively 
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FIG. 2: The crystal structure of ideal SrCrgGasOig. The light 
grey circles represent the oxygen atoms (we have restrained 
their number for clarity). The Sr atoms are not represented. 
The thick dashed lines show the typical hyperfine coupling 
paths of the gallium nuclei to various Cr sites through the 
oxygen ions. 



coupled to the Cr'^+ ions of the kagome bi-layer (Fig. |^). 
The present NMR study was carried out on different Cr- 
concentrations of SCGO (0.72 < p < 0.95). The compar- 
ative study of the samples shows that Ga(4/)-NMR can 
discern between the geometric frustration related prop- 
erties and the disorder related properties of SCGO: we 
are able to probe independently the kagome bi-layer sus- 
ceptibility and the effects related to dilution, and thus 
answer both questions on a firm experimental ground. 

The general outline of this paper is the following. Af- 
ter a detailed presentation of the samples used in this 
work, their characterization and a description of their 
magnetic couplings (Sec. |l^, we focus on the specificities 
of the Ga-NMR spectrum of SCGO: the gallium sites in 
SCGO, the nuclear parameters of the Ga(4/) resonance 
line and the optimization of its detection (Sec. [II), with 
a special emphasis on the new features brought by the 
study of the least diluted sample p — 0.95. We also show 
that NMR is a very refined technique for characterizing 
the amount of Ga/Cr substitution, especially in the low 
dilution limit. The following sections address the main 
points relevant to the frustrated physics. Compared to 
Ref. the study of the dilution effects is entirely new. 
After a brief overview of the raw NMR spectra and some 
experimental details (Sec. IV), we present in Sec. ^ the 
first part of our comparative study. By probing the sus- 
ceptibility on a local stand through the NMR shift, we 
identify the intrinsic, i.e. dilution-independent, suscep- 
tibility typical of the SCGO family — the archetype of 
geometrically frustrated compounds. We then discuss 
our experimental results in the light of the exis ting mod- 
els and calculations. The main topic of Sec. VI is the 
study of the dilution effects. Using the NMR width, we 
can track down in SCGO the effect associated with the 
dilution of the magnetic network and isolate its contri- 



bution to the macroscopic susceptibility. Moreover, we 
demonstrate that the impact of the site dilution is not an 
effect simply localized around the substitution site. As 
we discuss, this sets specific constraints on the theoretical 
models describing the dilution effects. We finally suggest 
an interpretation of our NMR results based on the gen- 
eral context of the AF correlated systems. A summary 



and the concluding remarks can be found in Sec. VII 



II. SAMPLES AND MAGNETIC COUPLINGS 

The Cr'^+ ions of SCGO occupy three distinct sites 
which are denoted 12k, 2a and 4/„i (Fig. |). The Cr(12fc) 
sites represent 2/3 of the total number of Cr'^+ ions and 
are arranged to form a kagome lattice. The two kagome 
planes are separated by a Cr(2a) triangular layer (1/9 of 
the Cr'^+ ions). The remaining 2/9 of the Cr'^^ ions oc- 
cupy the Cr(4/i,i) sites. A realistic description of the 
Cr-Cr magnetic couplings in SCGO showed that this 
compound has couplings similar to Cr203.L3 As was ev- 
idenced in Ref. p3|, SCGO is made up of only two mag- 
netic entities. IjThe kagome bi-layer, i.e. the kagome- 
Cr(2a)-kagome structure, with an average AF coupling 
of Jbi-iayer ~ 80 K as we establish in Sec. ^ 2) 
The Cr(4/ui)-Cr(4/i,i) spin pairs with an AF couphng 
of Jpair — 216(2) K, each pair being isolated from the 
others. The full structure is obtained by the stacking: 
spin pairs/kagome bi-layer/spin pairs/kagome bi-layer 
etc. The interaction between the kagome bi-layer and the 
spin pairs is small (~ 1 K). Since all the Cr(4/^,i) spin- 
pairs form non magnetic singlets at low-T (T <C Jpair, see 
Appendix), the low-T properties of SCGO are expected 
to reflect those of the kagome bi-layer network only. 

From the previous considerations, it is quite clear that 
the kagome bi-layer structure is more complex to model 
than a pure kagome system. However, SCGO is ideal 
in many regards. As mentioned in the introduction, the 
interactions in the magnetic hamiltonian other than a 
nearest neighbor interaction, are likely to modify the spin 
liquid nature of the ground state. These interactions are 
here extremely small compared to the Cr-Cr interaction 
in the kagome bi-layer. The anisotmpy of Cr''+ is « 0.08 
K (spins are therefore Heisenberg) ,Ea the dipolar interac- 
tion is ~ 0.1 K, and the next nearest neighbor interac- 



TABLE I: The characteristic susceptibility parameters for 
some of the SCGO samples studied. Me//(/^s) ^^'^ Qmacro are 
extracted by fitting the high-T behavior of Xm'acro — 150 
K) to a Curie- Weiss law oc [/Xgjj/(r -I- Qmacro)] 
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tions Jn.n.n SLie most likely | Jn.n.nl < 4 K, as in 

These microscopical details might not be totally uncorre- 

lated to the spin liquid-like behavior observed in SCGO. 

The present study was performed on a series of seven 
Cr-contents, corresponding to concentrations of p = 
0.72,0.81,0.89,0.90,0.91,0.93,0.95. All the samples are 
ceramics and were synthesized by a solid state reaction 
of SrCOg, CraOg and GaaOg in air at 1350° C, typi- 
cal for SCGO. The reaction products were checked by 
X-ray diffraction and by macroscopic susceptibility mea- 
surements. The physical parameters for Xmacro yielded 
results in agreement with literature (Tab. |). More re- 
fined characterizations were also performed. The p = 
0.81,0.89,0.95 samples were investigated by high resolu- 
tion neutron diffraction with the D2B spectrometer of the 
Insitut Laue-Langevinin in order to estimate the occupa- 
tion of the three Cr sites of SCGO, respectively pi2fc , P2a 
and Pif^- , otherwise non measurable by X-ray diffraction. 
The results are presented in Fig. The p-^ri-D.72 and 
0.90 samples were the object of /iSR studies.Eao Finally, 
the Cr-concentration of the samples was also checked di- 
rectly by Ga-NMR (detailed in Sec. 

The synthesis of samples with higher Cr-content than 
p = 0.95 failed. The X-ray diffraction showed a Cr203 
parasitical presence of > 0.5% (our detection thresh- 
old) in the final product oi p = 0.95 and all the addi- 
tional Cr203 introduced in the solid state reaction to 
reach higher concentrations than p = 0.95 simply did 
not react. This sets the limit of this synthesis method 
to the p = 0.95 concentration. Interestingly, this limit 
may be related to the fact that the Ga/Cr substitution 
in SCGO is not uniform, i.e. the Cr-concentration is 
not equal to p on all the sites. As shown in Fig. ||, 
the Ga/Cr substitution on the Cr(12A;) sites corresponds 
to stoichiometry {p « Pi2k), whereas in contrast the 




• Ref. 12, 33, 36 

□ This work 
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FIG. 3: The relative Cr-concentration on the 12k, 2a 
and 4/„i sites determined by the neutron measurements 
(pi2fc,P2a,P4/„i) versus the overall Cr-concentration (p) of our 
p = 0.81, 0.89, 0.95 samples (open squares) and of other sam- 
ples in the literature (closed circles). The Cr-concentration of 
the three sites follows stoichiometry when the symbols fall on 
the solid line. 
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FIG. 4: The Cr-Cr couplings as a function of the Cr-O-Cr 
angle (left) and the Cr-Cr distance (right) for: the a verag e 



Cr-Cr couphng of the SCGO kagop^ bi-lay(i*n(see Sec. VB), 
the SCGOi§r(4/„0-Cr(4/„0 pairsP CraOaP and a series of 
Cr-oxydes.EZI The stars stand for the Cr(12fc)-Cr(12fc) and the 
Cr(12fc)-Cr(2a) couplings of SCGO's kagome bi-layer which 
we determined in p-previous publication via a refined analysis 
of the NMR shift .Ef 



Cr(2a) sites are robust to substitution {p2a > p), the 
non-stoichiometric gallium occupying preferentially the 
Cr(4/^i) sites (p4/„i < p). We note that no gallium is 
present on the Cr(2a) sites of the p = 0.95 sample and 
that concomitantly the parasitical Cr203 problem occurs 
at this same Cr-concentration, an indication that the two 
phenomena may be related to the same chemical con- 
straint. 

A last point we would like to address, is the influence 
of the Ga/Cr substitution on the crystal parameters a 
and c. The largest structural change we observed is a 
w 0.004 A elongation of the Cr(12fc)-Cr(2a) distance in 
the p — 0.72 sample (2.975 A) compared to the p — 0.95 
sample (2.971 A). A study on a series of chromium-based 
oxides with Cr-O-Cr bonding angles ^ 90° (as in SCGO) 
showed that the Cr-Cr exchange constant decreases lin- 
early with increasing Cr-Cr distance (Fig. with a slope 
of AJ/Adcr-Cr « 450 K/A,EZl from which we extract a 
negligible variation of the SCGO exchange constants at 
most of A J « 3 K. We therefore expect the sizeable 
variation of the magnetic properties of SCGO with p to 
reflect only the effect of the lattice dilution. 



III. GA-NMR SPECTRUM 

The present NMR study was performed on the ^^Ga 
(69^ = 10.219 MHz/T, = 0.178 x lO'^* cm^) and 
the ^^Ga (^'^7 = 12.982 MHz/T, ^^Q = 0.112 x 10-^4 
cm^) nuclei of SCGO powder samples (7 and Q are re- 
spectively the gyromagnetic ratio and the quadrupolar 
moment of the nucleus). The gallium ions are present on 
two distinct crystallographic sites, which are designated 
by Ga(4/) and Ga(4e) (Fig. ^). A previous spectral anal- 
ysis carried out on the p = 0.90 sample in a sweep field 
set up at a radio-frequency of z^^/ — 131 MHz, success- 
fully assigned each site to the corresponding peak in the 



NMR spectrumBa There, it was shown that the Ga-NMR 
spectrum of both isotopes is actually the sum of three 
contributions: Ga(4/), Ga(4e) and an extra contribution 
related to the presence of non-stoichiometric gallium on 
the Cr sites, which we label by Ga(su6). 

As was pointed out in Ref. ^ the interest of Ga-NMR 
resides in the fact that the gallium nuclei are coupled to 
the neighboring magnetic Cr'^^" ions through a Ga-O-Cr 
hyperfine bridge (Fig. In particular, the ^^'''^Ga(4/) 
nuclei are exclusively coupled to the kagome bi-layer, 9 
from the 2 kagome adjacent Cr(12fc)-layers and 3 from 
the intermediate Cr(2a)-sites. Through Ga(4f)-NMR we 
are then able to probe locally the magnetic properties of 
the kagome bi-layer. The present NMR study is therefore 
devoted to Ga(4/)-NMR. 

Along with the hyperfine interaction mentioned here 
above, the nuclear hamiltonian of gallium in SCGO also 
bears an additional quadrupole interaction due to the 
electric field gradient (EFG) on the gallium sites. Fol- 
lowing the usual notations, the nuclear hamiltonian may 
be expressed as 



n = -h-/I ■{1+K)-H 
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[ill-I^+r^ill-lD] 



where H is the applied field, the principal axes of the 
magnetic shift tensor K are collinear with the direction 
of the nuclear spin operators /j,, ly, and I^, vq is the 
quadrupole frequency and < 77 < 1 is the quadrupole 
asymmetry parameter. The main focus of the spec- 
tral analysis presented in this section is to complete the 
study of Ref. ^ by determining the quadrupole contri- 
bution to the nuclear hamiltonian of Ga(4/), working 
in a frequency range (urf = 40.454 MHz) more appro- 
priate for the observation of the ^^'^^Ga(4/) resonance 
lines. We study here the spectrum of the p = 0.95 sam- 
ple which yields sharper quadrupole features than the 
p = 0.90 sample (or of any of the other samples studied). 
The knowledge of the quadrupole parameters reported 
m Tab. allows us, in the following section, to sepa- 
rate the magnetic and the quadrupole contributions to 
the Ga(4/) spectrum and, hence, to safely evaluate the 
magnetic properties of the kagome bi-layer. 

In the second part of this section we focus on Ga(sM6). 
We identify the contribution of Ga(su6) to the spectrum 



TABLE II: The quadrupole parameters of the 4e and 4/ gal- 
lium nuclei in SCGO. The ^^Ga isotope yields a stronger 
quadrupole interaction compared to the ^^Ga isotope, since 
the quadrupole moments are in a ratio of ^"Q/^^Q —^^ 
VQp^VQ = 1.589. 





^Vq (MHz) 


""^/Q (MHz) 


■n 


Ga(4e)" 


20.5(3) 


32.6(5) 


0.050(35) 


Ga(4/)'' 


2.9(2) 


4.6(2) 


0.005(6) 
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FIG. 5: Top panel. A field sweep of ^^Ga and ^^Ga at 
i^rf ~ 40.454 MHz. The site and transition assignments of 
each peak are described in the text. Middle panel. The satel- 
lite transitions in the ^^'^^Ga(4/) line. The satellite singu- 
larities of the 4e line are not resolved. Bottom panel. A line 
simulation of ^^•''^Gai4f) (solid line) and of ^^•''^Ga{4e) (dark 
grey area) performed on the basis of the quadrupole parame- 
ters of Tab. The simulated spectrum does not account for 
the Ga,{sub) contribution. 



acquired at Vrf — 40.454 MHz and, on the basis of this 
contribution, provide a new method to evaluate the Cr- 
concentration of the SCGO samples through Ga-NMR, 
especially well suited for the very low dilutions. 



A. Spectral analysis of Ga(4/) 

Neglecting the Ga-substituted sites, the Ga-NMR 
spectrum displays four sets of lines corresponding to the 
two isotopes distributed on both Ga(4e) and Ga(4/) 
sites. In a powder, for a given site and a given iso- 
tope, the lineshape results from the distribution of the 
angles between the field and the EFG principal axis. This 
yields singularities rather than well defined peaks — the 
so-called powder lineshape. 

The quadrupole interaction of gallium nuclei in SCGO 
is a consequence of the coupling of the nucleus to the 
EFG produced by the surrounding electronic charges, 
here mainly the oxygen ions. As shown in Fig. ^, Ga(4e) 
is surrounded by a bi-pyramid of 5 oxygen ions. This 
environment has no local cubic symmetry and results 
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into a strong quadrupole frequency as was evidenced 
in Ref. || (Tab. ||). In contrast with Ga(4e), Ga(4/) 
is surrounded by a nearly ideal tetrahedron of oxygens 
ions, only slightly elongated along the crystalline l^-axis. 
Since a regular tetrahedron has a local cubic symme- 
try, we expect the quadrupole frequency of Ga(4/) to 
be j^q(4/) ^ i^Q(4e). This agrees with Ref. pq where 
the Ga(4/) quadrupolar effects were not detected at a 
frequency of i>rf = 131 MHz. The asymmetry parameter 
rj, that quantifies the deviation of the EFG from axial 
symmetry, is in principle close to zero for both sites as 
they share a rotation axis along c?. 

We present in the top panel of Fig. |5| a field sweep 
spectrum, obtained at i/^/ = 40.454 MHz and T = 80 K 
for the p — 0.95 sample. The sweep, from 2.65 T to 4.80 
T, covers the spectrum of both isotopes and yields the 
expected features. We first focus on the ^^'^^Ga(4e) line. 
A powder sample simulation of the 4e-line is presented 
in the bottom panel of Fig. ^ to emphasize its contribu- 
tion to the spectrum. Each 4e-isotope, i.e. ^^Ga(4e) and 
^^Ga(4e), exhibits a resonance in a wide range of fields, 
with multiple quadrupole peaks and steps associated 
to the three nuclear Zeeman transitions 3/2 < — > 1/2, 
1/2 < — > -1/2 and -1/2 < — > -3/2. As evidenced in 
Fig. H, in this field window only the 1/2 < — > —1/2 res- 
onance transition is resolved yielding for each 4e-isotope 
two singularities known as the central line singularities 
(labeled respectively by CL— and CL+ in Fig. These 
two singularities are the boundaries of the central line 
splitting of width « z^q/27z/j./£j Compared to the work- 
ing frequency of i/rf = 131 MHz used in Ref. ^ the 
much lower frequency of i^rf = 40.454 MHz allows here 
to spread the 4e line over a wider field-range. The 
^^'^^Ga(4e) contribution to the spectrum is a nearly con- 
stant background of small amplitude compared to the 
^^Ga(4/) and ^^Ga(4/) lines — the two prominent peaks 
of Fig. — hence there is a better contrast between the 
sites. Furthermore, overlap between the Ga(4e) and the 
Ga(4/) lines is also minimized at this radio-frequency. 

We now focus on the ^^'''^Ga(4/) line of Fig. ^. In 
these favorable experimental conditions, it can be seen 
that both 4/ peaks are flanked by two shoulders, the so- 
called satellite singularities (labeled respectively by Sat— 
and Sat+ in the middle panel of Fig. p), associated with 
the 3/2 < — » 1/2 (Sat-) and the -1/2 < — > -3/2 {Sat+) 
nuclear transitions. The detection of the satellite singu- 
larities is a clear indication of the existence of quadrupole 
effects in the ^^^^'^Ga(4/) spectrum, implying also a cen- 
tral line splitting for Ga(4/). To emphasize this central 
line splitting, we present in Fig. ^ narrow field sweeps 
of both isotopes at T = 150 K. The sweep ranging from 
3.820 T to 3.960 T covers the ^^Ga(4/) line (bottom ab- 
scissa) whereas the sweep ranging from 3.026 T to 3.069 
T covers the ''-^Ga(4/) line (top abscissa). As shown, the 
two central lines match since we expanded the field win- 
dow of ^iGa(4/) by a factor ^^((9V7)/^^(0V7) = 3.209, 
i.e. by the ratio of the central line quadrupolar widths 
expected for the two isotopes. Fig. |6| also indicates that 



TABLE III: The field positions of the 4/ quadrupole singular- 
ities in Fig. ^ Hsing axe the fields where each singularity was 

the calculated fields where 

Tcalc 



found experimentally and Hlmg 

the singularities are expected. The missing Hsfn'g correspond 



to the field singularities used to perform the calculation. 



Sing. 


Using (T) 




71 rrcalc f'~r^\ 
^ sing \ / 


69 rrcalc ('~r^\ 
^ sing \ / 




3.025(2) 


3.838(1) 


3.026(2) 


3.838(2) 


4fCL+ 


3.036(1) 


3.865(1) 


3.035(2) 




^ jSat — 


2.924(9) 


3.628(8) 


2.921(9) 






3.137(7) 


4.066(8) 


3.138(7) 





the structure observed in the 4/ line is dominated, at 
these temperatures, by quadrupole effects. In particular, 
the anisotropy in the shift K, which could in principle 
yield additional structure to the line, is small, i.e. the 
shift is isotropic {K = K). This conclusion is further 
supported by a calculation, presented below, of all the 
field singularities. 

We now turn to the evaluation of the quadrupole pa- 
rameters of ^^''^^Ga(4/). This was carried out using the 
field positions Hsing oi the quadrupole singularities de- 
termined from Fig. and reported in Tab. [II which de- 
pend only on the three parameters z/q(4/), ^(4/) and 
provided the shift is isotropic: Hsing — H{i'Q,ri, K)SB 
We evaluated the 4/ quadrupole parameters from the 
singularity positions of 69Ga(4/)^"*-, ^^Ga.{Af)^^+ and 
69Ga(4/)^"*+. At T = 80 K, the shift is = 0.0285(5). 
The extracted quadrupole parameters are presented in 
Tab. H and are found in agreement with expectations. To 
confirm the isotropic nature of the shift, we performed a 
self-consistency test by calculating the field positions of 
all the quadrupole singularities in Fig. ^ This is read- 
ily done using the three parameters i'q(4/), ?7(4/) and K 
determined previously and the set of equations of Ref. Bfl. 
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FIG. 6: Field sweeps of the ^''Ga(4/) and the ■^^Ga(4/) lines 
at r = 150 K. The closed (open) symbols represent ^^Ga 
(^'^Ga). The top abscissa is scaled to the bottom by a fac- 
tor ^'^ {Q"^ h) {Q'^ h) = 3.209 to evidence the quadrupole 
structure in the lines. The hump at « 3.065 T is ^^Ga(sitfe). 
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As shown in Tab. [II, the calculated positions are in ex- 



cellent agreement with the experimental ones. 

Finally, concerning the quadrupole interaction in all 
the other SCGO samples studied, no appreciable change 
was evidenced in i/Q(4/) and »y(4/). In conclusion, the 
Ga(4/) central line is narrow enough in all the samples 
to allow to follow accurately the T-variation of the shift 
and of the linewidth. 



B. Gallium substituted on the Cr sites 

In Fig. we compare the ^^Ga spectra obtained for 
the p = 0.81, 0.89, 0.95 samples at T = 150 K and 
i^rf = 40.454 MHz. Care was taken to ensure the same 
experimental conditions. The field sweeps of Fig. ^ cover 
a narrow field window 6H (from 3.00 T to 3.15 T) cen- 
tered on the central line of ^iGa(4/) {H « 3.05T). 
Two features differ between the samples. 1) The line 
ai H 3.065T whose integrated area decreases with in- 
creasing p, which is therefore ^^Ga(sM&). 2) The con- 
stant background h{p) which decreases with increasing 



^^Ga(su6) and part of h{p) are therefore all the con- 
tribution p£ the non-stoichiomctric gallium to the ^^Ga 
spectruni.c3 

Since the non-stoichiometric gallium is present on the 
three Cr sites of SCGO, one would expect to detect three 
gallium lines and not a unique Ga(su6) line. To gain in- 
sight into this problem, we studied the spectrum of the 
parent and non magnetic compound SrGai20i9 {p — 0). 
The spectrum allows to evidence sizeable quadrupole ef- 
fects on the 12fc C'^uq = 5.1(1) MHz and = 0.47(3)) 
and on the Af^i ('^vq = 1.6(1) MHz and ij = 0.03(5)) 
substituted sites, otherwise-jaon measurable in the mag- 
netic compounds (p 7^ 0)r3 The quadrupole effects on 
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FIG. 7: A field sweep centered on the ^^Ga(4/) central line 
performed at T = 150 K. The solid line is a multiple gaus- 
sian fit. Inset, pnmr versus p. The open and closed symbols 
are respectively the Cr-concentration evaluated by direct in- 
tegration and by multiple gaussian fits. The dashed line is 
the expected value from neutron/X-ray difi'raction. 



the 2a site are negligible C^i^q ~ 0.3 MHz and 77 « 0). 
Turning back to the spectra of Fig. the ^^Ga(su6) 
line is the unresolved sum of the 2a line and of the 
12fc and 4/„i quadrupole central lines of substituted gal- 
lium, whereas the background b(p), in addition to the 
Ga( 4e) an d Ga(4/) satellites contributions mentioned in 
Sec. Ill A, contains the sum of the satellite contributions 



of the substituted gallium. 

In order to evaluate the Cr-concentration of our sam- 
ples using Ga-NMR {pnmr), we determine the intensi- 
ties lif and Isub{p) of respectively the ''^Ga(4/) central 
line and the ''^Ga(sM&) line of Fig. ^. To do so, we eval- 
uate the integrated area of both lines corrected for the 
transverse exponential relaxation (T2{sub) — 78(2)/is and 
^2(4/) — 50(l)^s). Isub{p) and 74/ are respectively pro- 
portional to the amount of gallium present on the three 
substitution sites and on the 4/ site. The lines are well 
separated in the p = 0.95 sample, so that Isub{p) and 
lif are evaluated separately. For the more diluted sam- 
ples, the intensity 14/ is fixed to the p — 0.95 value as 
the amount of gallium on the 4/ site is independent on 
Cr-concentration. Isub{p) was then evaluated either by 
direct integration of the spectrum over the field window 
6H and subsequent subtraction of both b{p)5H and 74/ , 
or by reproducing all the lines with a minimal 3 gaus- 
sian fit and an additional constant background b{p) (one 
gaussian for the ^^Ga(4/) line with an intensity fixed to 
the p — 0.95 value, and two gaussians for '^^Ga(sM6), the 
area of which yields I sub)- From the ratio Igub/Iif, we 
can extract pnmr- In the inset of Fig. ^ we show the vari- 
ation of Pnmr. with the concentration p determined by 
X-ray and neutron diffraction. Whatever the method em- 
ployed (integration or multiple gaussian fits), we notice a 
perfect agreement between the different characterization 
methods. Since pnmr is equal to p, we conclude that 
the actual chemical content corresponds to the nominal 
concentration. 

Incidentally, Ga-NMR proves to be a quite accurate 
tool for determining the Cr-concentration in SCGO, es- 
pecially for dilutions as low as 5% or less where the other 
techniques are beyond their sensitivity limit. For the 
samples with p > 0.95, we did not find any further de- 
crease of the Isub value which confirms locally that it is 
impossible to synthesize samples beyond the p = 0.95 
limit found through X-ray diffraction. 



IV. DILUTION AND T-DEPENDENCE OF THE 
GA(4/) SPECTRUM 

To appreciate the experimental evidence of the prop- 
erties we discuss in the next sections, we present here 
the dilution and the T-dependence of the raw ^^Ga(4/)- 
NMR spectra. The main relevant physical parameters of 
interest are the shift of the line, which as we establish al- 
lows us to probe directly the susceptibility of the kagome 
bi-layer, and the linewidth which in contrast yields infor- 
mation on the distribution of the internal local fields in 
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FIG. 8; Typical high-T field sweeps of the ''^Ga(4/) central 
line. 



the kagome bi-layer. We also briefly describe the analysis 
performed to extract these spectral parameters. 



A. Experimental data 

The Ga(4/) spectra were recorded in a wide T-range, 
from 5 K to 410 K. Up to 200 K, the spectra were ob- 
tained in a sweep field set up with a working frequency 
of = 40.454 MHz. For 150 K< T < 410 K, the spec- 
tra were recorded in a static field of « 7 T by a swent=. 
frequency variant of the field step spectroscopy method.c3 
Our data can be sorted in two distinct temperature do- 
mains, typically T > 100 K and T < 100 K. 

A typical series of high-T spectra (100 K< T < 180 
K) is presented in Fig. || for the p = 0.95 sample. The 
field window (from 3.00 T to 3.08 T) covers the central 
line of ^^Ga(4/). Fig. || is representative of the evolu- 
tion of Ga(4/) at T > 100 K for aU the samples stud- 
ied. The line shape is nearly constant with temperature, 
as dominated by T-independent quadrupole effects, and 
shifts towards the low fields with decreasing tempera- 
ture. The shift K, which is experimentally measured by 
the position of the ^^Ga(4/)'-^^+ singularity and perform- 
ing minor quadrupole corrections (see Sec. HI), therefore 
increases when the temperature is lowered. 

The T-dependence of both K and the width changes 
at low-T. We present in Fig. || spectra ranging from 5 K 
to 50 K for the p = 0.95 sample (top) and the p = 0.81 
sample (bottom). The field window covers the ^^Ga(4/) 
line, but is larger than in Fig. |, from 2.90 T to 3.15 T. 
Three important features should be noticed. 



1. The line now shifts toward the high fields and 
hence K is decreasing with the temperature. How- 
ever K does not reach the zero value, i.e. the 
^^Ga(4/) line is never centered at the reference 
value "^^Href = i^r//^^7 = 3.116 T. The compar- 
ison of the two series of spectra establishes that 



the change in the shift direction is not affected by 
the dilution. 

2. The width increases with decreasing temperature. 
The quadrupole structure that can still be noticed 
on the T = 50 K spectrum of the p = 0.95 sam- 
ple, is progressively washed away as temperature 
is lowered. At T = 20 K, the line is marked by a 
smooth symmetric broadening of gaussian nature. 
This broadening is more important for the p — 0.81 
sample than for the p — 0.95 sample, i.e. the width 
is sensitive to the dilution. The ^^Ga(su6) line is 
no longer resolved in this T-domain. 



3. 



A last feature we would like to underline is the 
abrupt decrease of the detected ^^'''^Ga(4/) nu- 
clear population at T < 15 K in both samples. 
The nuclear population, proportional to the inte- 
grated area of the 4/ line, actually decreases in 
the same way in all the samples studied (Fig. p^ . 
It therefore originates from a dilution-independent 
mechanism. The present result confirms the con- 
clusions of the preij-ious Ga(4/)-NMR study on the 
p = 0.90 sample.Lj There, the wipe-out of the in- 
tensity was assigned to originate from the intrinsic 
high dynamics of the kagome bi-layer spin system 
when T ^ K. This effect, expHcitly evidenced in 
SCGO by /iSR measurements,c3c3 is taken to be 
the signature of a spin liquid-like ground state. Be- 
cause of this abrupt decrease of the NMR signal. 
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FIG. 9: Low-T field sweeps of ^^Ga(4/) for the p=0.95 sam- 
ple (top panel) and the p=0.81 sample (bottom panel). The 
dashed line is a guide to the eye. 
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FIG. 10: The ®^''^iGa(4/) population detected by NMR. For 
each sample, the integrated area of the 4/ line is normalized 
by the integrated area measured in the 20 — 50 K temperature 
range, where it is found constant in temperature. The Ga(4e) 
contribution, constant in temperature, was subtracted. T2 
corrections were found quite small and do not affect estimates 
below 50 K. 



69^'^iGa(4/) is not a suitable probe for T < 10 - 15 
K. 

The comparison of the spectra of all the samples en- 
ables us in the next sections to easily separate the prop- 
erties related to the dilution from the intrinsic properties 
independent on the dilution. 



B. NMR line analysis 

Compared to Ref. we employed here a different 
method to evaluate the shift K and the width AH of the 
spectra of Fig. ^. (i) For all the samples, we choose a 
^^'^^Ga(4/) high-T reference spectrum dominated by T- 
independent quadrupole effects {Sq{H)) (ii) In order to 
recover the low-T broadening of gaussian nature, Sq{H) 
is convoluted by a normalized gaussian function G{H) 
of width AH. Since it is impossible to fit separately 
the substituted and non-substituted sites, we explicitly 
assume that for the highly diluted samples both sites 
broaden identically. This is likely the case in view of the 
hyperfine coupling paths. The low-T spectrum, 
is therefore simply reproduced by 



Sm{H) = a 



yo 



+ OC 



G{H-H')Sq{H' + Hs)dH' 



where Hs is the field-shift of Sq (H) , from which we ex- 
tract K. Since the ^^'^^Ga(4/) line lies on a constant 
background that results mainly from the contribution of 
the second gallium site Ga(4e) (see Fig. ||), the amplitude 
a and the constant background j/o are employed to read- 
just the relative spectral weights of the 4e and of the 4/ 



lines when the transverse relaxatioupli corrections varies 
differently on the two gallium sites.a We also restricted 
our fits to the left side of the line, very sharp at high-T 
and which has the major advantage of not being affected 
by the Ga/Cr substitutions. Concerning the values of K, 
we did not find any significant difference with the method 



employed in Ref. 32 



C. Magnetic contribution to the Ga(4/)-NMR line: 
shift and width 

We briefly recall here the relationship between the con- 
tribution of each gallium nucleus to the Ga(4/)-NMR line 
and its local magnetic environment in order to underline 
what can be exactly probed through Ga(4/)-NMR (more 
details can be found in Ref. p8| , p]j4^ ). 

Each Ga(4/) nucleus is coupled to its Cr(12fc) and 
Cr(2a) nearest neighbors (n.n.) through a Ga-O-Gr hy- 
perfine interaction with an hyperfine constant A.l3 We 
suppose that the susceptibility in the kagome bi-layer 
varies from Cr site to Cr site and label it, in a generic 
manner, by x- A Ga(4/) at site i will contribute to the 
NMR spectrum at a position depending upon the number 
of the n.n. occupied Cr sites and their susceptibility x- 
This corresponds to the shift i^T^*) in the NMR spectrum 
for a gallium at site i 



E 



Ax 



occupied n.n.Cr(12fc,2a) 

(the chemical shift is negligible). The average shift of the 
NMR line, that we label by K in the following sections, 
is simply related to the average susceptibility x over all 
the Cr sites. As we establish, the shift corresponds to 
the frustrated susceptibility Xfmstr^ so that K (x Xfmstr- 
Instead, the spatial distribution of K^^^ around K defines 
the magnetic width of the Ga-NMR spectrum and reflects 
the existence of a spatial distribution of x (a-s we detail in 
Sec. VI B , a spatial distribution of the hyperfine constant 



or a distribution related to the presence of Cr-vacancies 
in the nuclear environnement is negligible and cannot 
justify the low-T broadening observed). 



V. KAGOME BI-LAYER SUSCEPTIBILITY: 
EVIDENCE FOR CORRELATIONS EFFECTS 

As mentioned in the previous section, the shift K 
measures the average kagome bi-layer susceptibility that 
proves to be unobservable at low-T through the macro- 
scopic susceptibility measurements. We show here that 
the T-dependence of K gives a sharp evidence for a maxi- 
mum in the susceptibility Xfmstr around 50 K, whatever 
the amount of dilution. The comparison of our results on 
the SCGO compound with the different models existing 
in the literature, does not allow us to totally validate or 
invalidate them. Our observations rather favor an image. 
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FIG. 11: Top panel. K versus T down to 15 K. A minor 
second-order quadrupole correction has been performed. In- 
set, versus T for p — 0.95. Bottom panel. K versus T 
plotted using reduced units. Ccw and Qnmr are respectively 
the Curie- Weiss constant and temperature extracted from the 
high-T Curie- Weiss fit to described in the text. 



supported by recent neutron diffraction data, where the 
short range magnetic correlations play a central role. 



A. An intermediate temperature scale 



in ii' a well known property of the susceptibility of the 
geometrically frustrated compounds: the Curie- Weiss be- 
havior continues to subsist at temperatures T <C Qnmb.- 
Assuming the Curie- Weiss value of the kagome bi-layer 
to be (z) Jbi-iayerS{S + l)/3kB, whcrc (z) = 5.14 is the 
average number of nearest neighbors for a Cr site of the 
kagome bi-layer as the Cr-environments of the Cr(2a) 
and the Cr(12fc) sites differ, we find an average exchange 
constant of Ju-iayer ~ 80 K. As shown in Fig. ^, this 
value is close to the exchange constant J = 76(3) K ob- 
served in Cr203 where the Cr'^+ ions have a local octa- 
hedral environment similar to Cr(12A:).E3 We also gained 
more insight on the couplings of the kagome bi-layer 
through a mean- field analysis of the high-T shift, detailed 
in Ref. ^ that allowed us to evaluate the Cr(12fc)-Cr(2a) 
and the Cr(12A;)-Cr(12fc) couplings (their values are re- 
ported on Fig. Turning back to the NMR Curie- 
Weiss temperature, Qnmr is of the same order as the 
macroscopic Curie- Weiss temperature, although smaller 
with an increasing difference at higher Cr-concentrations. 
This difference is related to the fact that the susceptibil- 
ity of the Cr(4/i,i)-Cr(4/t,i) isolated spin pairs also con- 
tributes to XmacroJ but is abscut in K since the Ga(4/) 
nuclei probe only the susceptibility of the kagome bi-layer 
(Sec. [VB] ). 

The linear behavior of if for T ^ 8 nmr and the de- 
viation from this behavior that results in a maximum 
in K versus temperature, is a common feature of the 
three samples. It then stems from a physics robust on 
dilution, therefore related to an intrinsic property of 
the kagome bi-layer. Tmax is a new temperature scale 
for SCGO since it differs substantially from the two 
known characteristic temperatures, the freezing temper- 
ature Tg {Tmax/Tg « 10) and the Curie- Weiss tempera- 
ture {Tmax/QNMR « 0.1). 



The T-dependence of K for the p = 0.81,0.89,0.95 
samples is presented in Fig. The figure is the quanti- 
tative evaluation of the line shift of Fig. || and Fig. ^. K 
increases following a Curie- Weiss law (see below) up to a 
temperature of Tmax = 40— 50 K, where it reaches a max- 
imum. Below Tmax, K decreases, but does not reach for 
the observed temperatures (T > 15 K). The decrease in 
K is slightly more pronounced with increasing dilution. 
For the p = 0.81 sample K{T = 15 K)/K{Tmax) « 0.86, 
instead of 0.92 for the p = 0.95 sample. Although for 
T < 15 K the measurement of K is less significative 
since a wipe out of the intensity occurs, we observed 
in the p — 0.90 sample only a decrease of K{T = 7 
K)/K{Tmax) ~ 0.77.EI 

The typical T-variation of K^^ is presented in the 
inset of Fig. Ill] for the p = 0.95 sample. The high- 
T variation of R-^ (T > 100 K) is linear for aU the 
samples, suggesting a Curie- Weiss behavior. A linear 
extrapolation to = yields the Curie- Weiss tem- 

perature as determined by NMR (Qnmr)- The values 
of Qnmr are 453(30) K, 469(27) K and 484(25) K for 
the p = 0.81, 0.89, 0.95 sample respectively. We find here 



B. Evidence for a two component macroscopic 
susceptibility 



We compare in Fig. |lj the T-dependence of K (left 
y-axis) and of Xmacro (right y-axis) for the p = 0.95 sam- 
ple. The discrepancy between the low-T behavior of K 
and Xmacro Quitc clear since instead of a maximum 
in the susceptibility, Xmacro exhibits a Curie-like law at 
low-T. The difference between K and Xmacro is also per- 
ceptible directly on the spectra of Fig. ^ where an arrow 
indicates the approximate position where the maximum 
of the ''^Ga(4/) line should occur if cx Xmacro- In- 
stead of a model where Xmacro has a single component, 
our results establish that Xmacro yields at least two dis- 
tinct contributions. The first one is Xfmstri spatially 
uniform over the kagome bi-layer and that is reflected in 
K (Sec. W^ . A second one, non-uniform over the Cr 
sites, that is necessary to explain the Curie-like upturn 
observed at low-T. We label this contribution by Xdef^ 
since m Sec.0 we establish that Xdef is the susceptibility 
of the magnetic defects generated by the dilution of the 
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FIG. 12: K (left) and Xmacro {right) for the p=0.95 sample. 
K and Xmacro do not follow exactly the same law at high-T, 



since Xmacro piobes also the susceptibility of the Cr(4/^i)- 



TABLE IV: The fitting parameters ^^^.^.(4/„i), Cdef and A 
for the n — 0.81, 0.89, 0.95 samples used to reproduce Xmacro 
by Eq.(0). 



p 


tJ-cffWvi){lJ-B) 


Cdcf (emu K/mol) 


A X 10^ (Oe/^ig) 


0.81 


4.0(6) 


0.51(12) 


3.4(4) 


0.89 


4.5(4) 


0.45(10) 


3.8(4) 


0.95 


4.3(2) 


0.25(5) 


3.9(2) 



Tab. where the C^e/ value is dominantly determined 
by the very low-T region fit of Xmacro (Sec. VI). Finally, 
A is nearly constant with p, confirming the assumption 
of Sec. H that the site dilution does not alter in a signifi- 
cant way the couplings in SCGO, hence that the Ga-O-Cr 
hyperfine interaction is also constant with the dilution. 



Cr(4/t,i) spin pairs. 



Discussion 



kagome bi-layer. The comparison of Fig. |lj establishes 
on an experimental ground the conjecture by Schiffer et 
al. of a two component macroscopic susceptibility^ 

The macroscopic susceptibility probes both contribu- 
tions at low-T. If we also take into account the suscepti- 
biHty of the Cr(4/t,i)-Cr(4/i,i) isolated spin pairs {Xpair)i 
the three contributions to Xmacro ^'^'^ 

Xmacro Xfrustr ~^ Xpair ~^ ^def ("^) 

The dominant Curie susceptibility of the defects Xdef 
at low-T prevents to probe Xfmstr through macroscopic 
measurements. 

Using Eq.(0), we can fit the T-dependence of Xmacro 
given some simple remarks and minor assumptions. The 
T-dependence of Xfmstr is known through K, and for 
this reason the fit was performed in the 15 K< T < 350 
K range, where we have a full intensity in the NMR sig- 
nal. This T-range of the fit leads to an inaccuracy in the 
determination of a possible low-T Curie- Weiss tempera- 
ture of Xmacroi hcncc wc assumc a pure Curie component 
Xdef = CdcfjT to fit the low-T upturn. The suscep- 
tibility Xpair is quantitatively derived analytically (see 
Appendix) . The only unknown parameters of Eq. (n|) are 
the effective moment ^g^j(4/t,i) of the Cr'^+ ions olXpain 
the Curie constant Cdcf of the defects and the hyperfine 
constant A. 

The fit to Xmacro ^hc -p = 0.95 sample is presented 
in Fig. ^ and reproduces correctly the evolution of the 
experimental macroscopic susceptibility in the T > 15 K 
range. In Tab. IV we have reported the fitting parame- 
ters extracted from the three samples -p — 0.81, 0.89, 0.95. 
The effective moment is close to the 3.87/1^ value ex- 
pected for a Cr^+ ion, and the Curie contribution Cdcf is 
found to decrease with increasing -p. The values of Cdcf 
presented here, even if extracted from a fit which does 
not cover the 5 K< T < 15 K range where Xdcf domi- 
nates, nevertheless are in good agreement with those of 



We compare here our experimental data for the intrin- 
sic susceptibility x fmstr '^'^ frustrated lattice reflected 
in K to the calculations performed for the kagome and 
the pyrochlore networks through different approaches. 

Within a classical frame where the minimization of the 
Heisenberg hamiltonian is constructed on the basis of 
non magnetic triangles and tetrahedras, different authors 
have simulated the susceptibility of the kagome 
pyrochlore network by Monte-Carlo calculations 
All these models suggest a Curie- Weiss susceptibility 
which extends to low-T (T <C O), in agreement with the 
high-T behavior of if (T > 100 K). Although these mod- 
els do not predict a maximum in the susceptibility versus 
temperature, all agree with our results on the absence of 
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FIG. 13: The three contributions to Xmacro obtained through 
the fit described in the text: the kagome bi-layer suscepti- 
bility (x/rustr)i the Cr(4/„i)-Cr(4/„i) spin pairs susceptibil- 
ity (Xpair) s-iid the Curie susceptibility induced by dilution 



(X 



dcf 



C del IT). 



12 



a Curie upturn at low-T as is detected in Xmacro- 

The susceptibility for the S = 1/2 kagome network 
has been calculated by aa exact diagonalization of the 
Heisenberg hamiltonian.O As mentioned in the intro- 
duction, the ground state is here constructed on non 
magnetic spin singlets and a "gap" (A) is predicted be- 
tween the singlet ground state and the magnetic spin 
triplet state. The calculated susceptibility hence varies 
as ~ exp{^A/T) at low-T, with a maximum in the sus- 
ceptibility occurring at A w O.IJ. However, it seems dif- 
ficult to reconcile this maximum with our results. First, 
we would expect to observe in K a much sharper de- 
crease in temperature than what is observed experimen- 
tally. Second, the energy scale of A itself does not agree 
with our data. In fact, taking J « 80 K, the maximum 
of K falls at T^ax ~ 0.5 J > A (note that A should he. 
weaker than O.IJ for S = 3/2 spins or even disappear) .ca 
The interpretation of the observed maximum for SCGO 
has to be researched elsewhere than in the gap predicted 
by this model. We stress that the Ga(4f)-NMR intensity 
loss does not allow us to probe deeply the predictions of 
this model which places the gap in the T < 15 K range. 
But, clearly, our data reveal that other parameters have 
to be taken into account for T 3> A to explain the origin 
of the maximum in if vs T at ~ 0.10. 

In this context, it appears natural to seek an inter- 
pretation in more conventional terms. The maximum in 
K versus temperature indicates an increasing magnetic 
rigidity of the kagome bi-layer spin network. As in all 
the d = 2 AF systems, we suggest that the maximum is 
an experimental signature of a small reinforcement of the 
magnetic correlations. The quasi-constancy of K for all p 
below this maximum (only a slight decrease is observed) 
is therefore an indication that the correlations are short 
ranged and concomitantly do not depend on the dilution. 

Let us recall the typical characteristics evidenced in 
the susceptibility of d = 2 AF networks with Heisenberg 
spins as reviewed in Ref. At high-T, the d = 2 sus- 
ceptibility follows a Curie- Weiss law, but exhibits a max- 
imum in the susceptibility at a temperature T^ax ~ ©, 
which is a consequence of the development of magnetic 
correlations. However, in contrast with d = 3 systems 
where a long range order is established at Tc ~ O, a 
d — 2 system has a finite correlation length at T^ax ~ © 
that eventually diverges when T ^ K. The low dimen- 
sionality of the lattice inhibits the full development of the 
critical fluctuations. In the case of the kagome bi-layer 
of SCGO, we suggest that the geometric frustration is 
even more efhcient in preventing correlations to develop, 
moving the maximum in Xfrustr down to T^ax ~ 0.18. 

On the experimental side, new neutron diffraction mea- 
surements further support this point of view. As men- 
tioned in Sec. I, neutron studies on SCGO have estab- 
lished that the magnetic correlations in the SG phase 
(T < Tg) are short ranged. C. Mondelli et al. observed 
more recently the neutron diffraction spectrum in a wider 
temperature range— than the previous experiments (1.5 
K< T < 200 K)E] As shown in Ref. 1^, the charac- 



teristic diffuse peak typical for short range correlations 
(^ « 2 X dcr-Cr) develops at T < 60 K, precisely in the 
same temperature range where the maximum in x fmstr 
occurs. This maximum is therefore the signature in the 
susceptibility of the development of these short range cor- 
relations. The feeble decrease in Xfrustr observe sim- 
ply reflects the finite value of ^ at T < Tmax ■ 

A recent quantum mean field theory study also points 
at this direction.Ed The basic idea here is to explicitly ac- 
count for the magnetic correlations by constructing the 
ground state not on interacting spins, but on interact- 
ing triangles or tetrahedras of spins. Within this frame- 
work, Garcfa-Adeva et al. calculated the susceptibility 
of the kagome and the pyrochlore networks for various 
spin values, providing a reference for comparisons to our 
Ga-NMR shift K (their mean-field theory does not ac- 
count for the existence of a spatially non uniform suscep- 
tibility responsible for the NMR-.linewidth broadening, 
see Sec. VI). In a recent paper J£j they also studied the 



influence of the lattice dilution on the susceptibility of 
these networks. This theory seems to capture most of 
the main features of our data: the calculated suscepti- 
bility for the diluted kagome and pyrochlore networks 
exhibits a maximum at T,nax ^ 0.15O followed by a fee- 
ble decrease of the susceptibility in the 0.056 — 0.159 
temperature range. Although the order of magnitudes 
found are very encouraging, as we underline in the bot- 
tom panel of Fig. |l^, more experimental and theoreti- 
cal studies are needed here to completely validate this 
model. Indeed, Garcfa-Adeva et al. predict that the 
susceptibility's maximum should be progressively washed 
out with increasing dilution and eventually disappear (in 
a S" = 5/2 kagome lattice for a ^ 10 — 15% dilution and 
in a = 3/2 pyrochlore lattice for a ~ 20 — 30% dilu- 
tion), a feature that somewhat contradicts our observa- 
tions. Furthermore, as evidenced in the bottom panel of 
Fig. |ll|, the T-dependence of K for T > Tmax is dilution- 
independent. This would only agree with the suscepti- 
bility calculated for a diluted pyrochlore lattice, but for 
T < Tmax the more pronounced decrease in K we observe 
is opposite to what is predicted for the pyrochlore lat- 
tice. Un fortu nately the Ga-NMR intensity loss at T < 15 
K (Sec. IV A ) does not enable us to further verify their 
model, in particular to verify whether there is a Curie- 
like upturn in the susceptibility for T <C 0.050. We con- 
clude then that for an adequate comparison between data 
and theory, the influence of the kagome bi-layer structure 
should be taken into account in the calculations. Also, 
the non random distribution of the non magnetic vacan- 
cies in SCGO's kagome bi-layer might induce quantita- 
tive differences between the calculated and experimental 
susceptibilities, as pointed out in Ref. ^ 

The general picture drawn by our overview of the ex- 
perimental and theoretical data is that even short range 
spin-spin correlations clearly play a major role in SCGO, 
or more generally in all geometrically frustrated com- 
pounds. Future studies are required in this direction to 
further uncover the low-T properties of these systems. 
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TABLE V: The parameters of the two component fit to Xmacro described in the text. At low dilution, the fit yields a stronger 
error. 



_P_ C (emu K/mol) 9 (K) Cd,f (emu K/mol) 9^^/ (K) Cj^f/C (%) 

0.72 14.0(2) 467(19) 0.53(3) 2.3(1.0) 3.8(2) 

0.81 17.6(3) 650(15) 0.39(3) 0.9(1.0) 2.2(2) 

0.89 20.5(4) 752(17) 0.34(3) 0.6(1.0) 1.7(2) 

0.90 22.0(6) 792(25) 0.24(3) 0.5(2.0) 1.0(2) 

0.95 23.5(7) 826(22) 0.15(3) 0.2(2.0) 0.6(2) 



VI. DILUTION EFFECTS 

In this section, we focus on the dilution-dependent sus- 
ceptibihty revealed by the NMR lincwidth at low-T. We 
show that the low-T paramagnetic behavior observed in 
Xmacro Originates from the dilution of the kagome bi- 
layer. This feature is observed in nearly all the geometri- 
cally frustrated AF compounds anjl-appears as a generic 
behavior of this class of materials .Eil 

Naturally, the macroscopic susceptibility could allow, 
in some way, to appreciate the dilution effects through its 
low-T behavior. Though, the multiplicity of the Cr sites 
in SCGO complicates the analysis. The Ga(4/)-NMR 
width not only enables to evidence that the vacancy of a 
spin on the network, i.e. the dilution, generates a param- 
agnetic defect, but also allows to shed light on the more 
fundamental question concerning the nature of the de- 
fect. 



A. Low-T macroscopic susceptibility 

In Fig. |lj we present the low-T macroscopic suscepti- 
bility of most of the samples studied. The low-T para- 
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FIG. 14: The low-T macroscopic susceptibility of some of the 
SCGO samples studied. The solid line is the two component 
fit (Eq.(|)) to the p = 0.90 sample. 



magnetic upturn of Xmacro observed at all the concen- 
trations, increasing in an appreciable way with growing 
dilution. This establishes that the origin of this contri- 
bution is related to the Ga/Cr substitution. The non 
magnetic gallium vacancy on the network must there- 
fore induce a perturbation which affects the neighboring 
magnetic Cr'^"'' ions. A paramagnetic defect is then gen- 
erated, whose susceptibility Xdef is reflected in the low-T 

behavior of Xmacro- 

In order to determine the exact low-T dependence of 
the susceptibility associated with the defects, the con- 
tribution to Xmacro of both Xfrustr and Xpair should in 
principle not be neglected. Since Xfrustr cannot be deter- 
mined through NMR below 15 K, we follow Ref. ^ and 
conveniently fit Xmacro by a two component expression 



C 



Cd. 



'ef 



T+e T+Q 



def 



(2) 



Anyhow, the corrections from Xfrustr and Xpair are small 
at low-T compared to Xdef and do not affect significantly 
the analysis. The first Curie- Weiss term roughly takes 
into account the contribution from Xfrustr and Xpair ^ 
which dominate at high-T. The second term, the more 
relevant for this section, also Curie- Weiss, quantifies the 
contribution of Xdef- 

The fits, shown by a solid line in Fig. ^for the p = 0.90 
sample, were performed in the 5 K < T < 350 K range 
for all the samples and reproduce correctly the low-T be- 
havior of Xmacro- The fitting parameters are presented in 
Tab. and are in agreement with Ref. 3ll The values of 
C and Q yield an overestimate of ~ 10 — 20% compared 
to those of Tab. |. The values of Cdef are in agreement 
with the analysis of Sec. 0. As expected, Cdef decreases 
with the dilution and most importantly is found to be 
only a weak fraction of the constant C, Cdef/C ^ 1% 
(Tab 0). Finally, the Curie- Weiss temperature O^e/, 
which is an indication of the average interaction between 
the defects, is negligibly small within error bars except 
for the p = 0.72 sample for which it is at most 2 K. 

Strictly speaking, it should be finally noticed that 
Xmacro not Sufficient on its own to establish experi- 
mentally whether the defects stem from the dilution of 
the kagome bi-layer. Indeed, the Ga/Cr substitution on 
a Cr(4/^i) site can in principle break a spin pair and 
free a paramagnetic spin, that can then contribute to 
the paramagnetic upturn of Xmacro- In contrast with 
Xmacro^ Ga(4/)-NMR is not sensitive to Cr(4/i,i) and, as 
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presented in the next subsection, allows a better under- 
standing of the dilution effects in SCGO. 



B. Dilution effects through the Ga(4f)-NMR width 

In Fig. |l5| we present the T-dependence (10 K< T < 
60 K) of the low-T width of ^^''^^Ga.{Af) {^^''^^AH) for 
five of the seven samples studied. We recover here the 
properties mentioned about the raw spectra: ^^'^^ATJ 
increases as the temperature drops and is very sensitive 
to the dilution, in contrast with K. The low-T behavior 
of ^^''^^AH bears strong similarities to the one of Xmacro- 
The perfect scaling of the two isotopes widths ^^ '^^AH 
normalized by the reference field ^^'^"'^iJr-e/ '^rf /^^''''^'y) 
— explicitly shown for the p = 0.95 sample in Fig.^ — 
underlines the magnetic origin of the low-T broadening 
(^69,71^^ 0^69,71 Href), p, agreement with the prior study 
on the p = 0.90 samplcEj 

The fact that the ^^'''^Ga(4/) linewidth increases 
rapidly at low-T and with the dilution, whereas the shift 
varies little, establishes the existence of a susceptibility 
that is spatially inhomogeneous over the kagome bi-layer, 
due to the defects associated with the Ga/Cr substitu- 
tion. Since the Ga(4/) nuclei are only coupled to the 
kagome bi-layer, the defects probed by Ga-NMR are nec- 
essarily localized in the kagome bi-layer. 

The effect associated to the dilution is not reduced to 
a mere suppression of a Cr site without any magnetic sig- 
nature over the neighboring spins (i.e. a magnetic defect) 
as indicated by the large modification of the ^^'^^Ga(4/) 
lineshape at low-T. Indeed, in the opposite case, one 
would expect a lineshape which results from the gal- 
lium nuclear spins coupled to various Cr-environments 
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FIG. 15: AH versus r for some samples studied. ^^'''^AH 
is extracted from the V rj = 40.454 MHz spectra following 
the fit described in Sec 
reference field '^^'"H,, 
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FIG. 16: Top panel: Cdef extracted from Xmacro iWt) and 
Cdef/NMR [right) extracted from the NMR width ^^'"^^ AH 
versus dilution (1 — p). The dashed line is a guide to the eye. 
Bottom panel: The SG temperature Tg of our SCGO samples 
versus dilution. 



all bearing the same susceptibility. One would therefore 
expect the linewidth to scale with the susceptibility, then 
_ftr, which varies only by 20% from 60 K to 15 K, whereas 
the linewidth increases by one order of magnitude. For 
this very same reason, a spatial distribution of the hy- 
perfine constant, which yields also a width cx K , cannot 
justify the broadening observed and also has to be ruled 
out. 

The T-dependence of ^'^ ''^AH is correctly reproduced 
by a Curie- Weiss law of the form C^^f /nmrI{T + 
Qdef/NMR), where Cdef/NMR is a constant and 
® def /NMR a Curie- Weiss temperature. The fits are pre- 
sented in Fig. ^as solid hues. The values of Qdef/NMR 
are extremely low, at maximum of w 4 K and vary ran- 
domly with dilution. Although we cannot exclude com- 
pletely a Curie- Weiss behavior, especially at high dilu- 
tion, we may consider the evolution in temperature of 
69,71^^ as very close to a Curie law. The values of 
Cdef/NMR are plotted as a function of (1 — p) in the 
top panel of Fig. |l^ (right y-axis) along with the valuesj 
of Cdef extracted from Xmacro by Eq- @ (left y-axis) £3 
As shown, the dependencies of both Cdef and Cdef /nmr 
on (1 — p) perfectly match, suggesting that ^^ "^^AH and 
Xmacro probc the same susceptibility at low-T, i.e. Xdef- 
This is further underlined in Fig. |l^ by the linear re- 
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FIG. 17: ^^''^'^AH versus Xmacro some of the samples stud- 
ied (T is an implicit parameter). 



lationship between ^^''^^ AH /Href and Xmacro at low-T. 
Since the defects probed by Ga-NMR are located in the 
kagome bi-layer, Fig. ^ and Fig. |l^ indicate that the 
defects on the Cr(4/i,i) sites yield a contribution scaling 
with the one from the kagome bi-layer or little contribute 
to Xmacro- 1^ vicw of the quantitative analysis of Cdef 
presented in Rcf. ^ the latter explanation is the more 
plausible. Finally, we note in Fig. |l^ that Cdef/NMR and 
Cdef are linear with (1 — p) at low dilutions {p > 0.90), 
and become progressively sub- linear with (1— p) at higher 
dilutions. This deviation is directly perceptible on the 
spectra of Fig. H, where the FWHM of the p = 0.81 sam- 
ple is clearly not « 4 times that of the p = 0.95 sample. 



C. Discussion on the nature of the defect 

1. Models for the localized defects 

To our knowledge, two models describe the dilution ef- 
fects in a geometrically frustrated system. Both are con- 
structed in a classical approach where the ground state 
corresponds to a minimization of the exchange energies 
on triangles (or tetrahedras) of spins. 

IChe most recent model was proposed by Moessner et 
a/.l£3 The basic idea here is that a paramagnetic mo- 
ment in the kagome network is generated when two va- 
cancies are adjacent on a triangle (3 adjacent vacan- 
cies are needed for the pyrochlore network). The de- 
fect corresponds in this case to a unique spin of the net- 
work, and at low dilution the defect's susceptibility is 
Xdef ~ (1 for the kagome network (~ {l—p)-^/T 

for the pyrochlore network). This prediction is in strik- 
ing contrast with the linear variation of the NMR width 
on dilution reported in Fig. |l^. To reconcile our data to 



the theory of Ref. |4^, it was suggested in Ref. ^ that 
the magnetic broadening of the width may reflect a dis- 
tribution of random local fields (cx Xdef)^ due to the cou- 
pling of each gallium nucleus to the 12 neighboring Cr 
sites of the kagome bi-layer. Basically, the NMR width 
Ai7 would then probe the mean deviation of the sus- 
ceptibility on these 12 sites (and not simply Xdef )i thus 
IS.H (X Y^(l — pY/T. Although appealing, this scenario 
is quite unlikely since the macroscopic susceptibility mea- 
surements which are insensitive to the effect of a local 
summation of fields, do not yield a quadratic (1 — p)"^ 
variation on dilution. 

The second model — purely empirical — was proposed 
in Ref. There, a defect is generated by a vacancy 
which perturbs the neighboring spins, uncorrelating them 
from the rest of the spin network. These uncorrelated 
spins were baptized "orphan" spins. Two distinct spin 
populations then coexist: the orphan spins with a para- 
magnetic behavior making up the defect, and the frus- 
trated spins outside the defect with a susceptibility equal 
to Xfrustr- In order to examine the implications of the 
orphan model on the Ga(4/)-NMR low-T spectrum, we 
assume, as a starting point, that the defect is built only 
on the orphan spins nearest neighbors (n.n.) of the va- 
cancy (Fig. |l^). We consider the important case of the 
p — 0.95 sample where the dilution is weak enough to 
allow refined conclusions and where the Ga/Cr substitu- 
tion occurs only on the kagome layers (the Cr(2a) site is 
not substituted, see Fig. ^. A simulation allows us to 
model the expected lineshape. To do so, we start by ran- 
domly suppressing 5% of the spins on the kagome network 
to account for the existence of non magnetic vacancies. 
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FIG. 18: The simulated sweep field central line of '^^Ga(4/) at 
T = 20 K and T = 10 K {i^rf = 40A5'iMHz) if orphan spins 
n.n. to the vacancy were present in the kagome bi-layer of 
the p = 0.95 sample. Also shown is the experimental central 
line of ^^Ga(4/) recorded for p = 0.95 at 20 K (open circles). 
The area of the experimental and of the simulated lines are 
normalized. 
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The n.n. spins of the vacancies are given a suscepti- 
bihty oc j[def) /T where iJ,^fj{def) « l-S/i^ (evalu- 
ated later in the section) to account for the orphan spins, 
whereas the remaining Cr sites a susceptibility X f rustr- 
Next, we classify the various nuclear populations along 
their Cr-environments. To construct the spectrum, each 
population is associated to a quadrupole line simulated 
with the parameters of Tab. with a shift reflecting 
their Cr-environment and an intensity weighted by the 
population size. Finally, the Ga(4f)-NMR simulated line 
is obtained by summing the spectra associated to all the 
nuclear populations. The line simulation for T = 20 K 
and T = 10 K is presented in Fig. |l^. As shown, neither 
the broadening, nor the simulated lineshape agree with 
our ^^Ga(4/) low-T experimental spectrum. Indeed, in 
the p — 0.95 sample, « 35% of the gallium nuclei do not 
probe an orphan spin, hence « 35% of the spectral weight 
yields a shift oc Xf rustr (short dashed line in Fig. [18| ). 
The majority of the gallium nuclei probe at least one or- 
phan spin, so that « 65% of the spectral weight yields a 
shift dominated by a 1/T contribution at low-T (dashed 
line in Fig. p^ . Independently on the value of fJ-^fj, one 
would therefore expect to observe « 65% of the spec- 
tral weight to shift toward the low fields with decreasing 
temperature whereas the remaining part is little shifted. 
This would yield a non-observed asymmetric broadening 
of the NMR line with decreasing temperature. 

Clearly, a model built on orphan spins n.n. of the va- 
cancy is unsatisfactory and the perturbation generated 
by the vacancy must be extended in space. Can we then 
imagine the defect to be made up of orphan spins which 
are first, second, third etc. neighbors of the vacancy? 
Since the spin-spi&jnteraction is AF this situation is very 
unlikely to occur .t3 As we show below, an extended stag- 
gered perturbation constructed on an AF interaction is 
conceivable. 



2. A defect built on AF correlations 

In Sec. we showed via the observed NMR shift that 
Xf rustr exhibits a maximum in temperature, a common 
feature for the AF correlated systems. It is therefore nat- 
ural to seek an interpretation for the origin of the para- 
magnetic defect also in the general context of the AF 
correlated systems. Indeed in the AF systems such 
the d — 1 spin chainsHLthe quasi— c? = 2 spin ladders,Ea 
and the d = 2 cupratesj£3 it is now well established that a 
vacancy (or a magnetic impurity) generates a long range 
oscillating magnetic perturbation and creates a paramag- 
netic component in the macroscopic susceptibility. The 
symmetric broadening of the NMR line observed in these 
systems is to be related to the oscillating character of 
the perturbation. We propose that the dilution effects of 
the SCGO kagome bi-layer can be described by the same 
physics of these correlated systems. 

We recall first some basics concerning this model. The 
presence of a vacancy in the AF correlated network of 



spins develops a magnetic perturbation, which, in a gen- 
eral fashion, affects a spin at the lattice position 'r by 

m(r) x'(r)Href 

where m(r), x'(^) ^^'^ Href are the magnetization, the 
susceptibility and the applied magnetic field (constant 
over r). For an AF system, x'(^) is peaked at a given 
vector Q. As an example, we choose x'{^) to have a 
gaussian shape 

X'{r) = x*g{r ■ Q)exp(-rV4f ) 

where the function g(r ■ Q) is oscillating and periodic 

in y with a periodicity related to Q. ^ is the spin-spin 
correlation length and x* the amplitude of x'(r'). This 
translates into an oscillating polarization of the spin net- 
work with a periodicity of ~ and damped over ~ ^. 

In Fig. |l^ we illustrate this with a polarization cal- 
culated for the kagome network. In the absence of low 
diluted single crystals necessary to determine the exact 
form of x'(^) and the nature of the vector Q through neu- 
tron diffraction experiments, we assume a gaussian shape 
for x'(^) and take the reciprocal vector Q — (1/3, 1/3) 
of the so-called kagome v^x a/3 spin configuration. By 
selecting this vector, we partially account for the effects 
related to frustration. 

The staggered shape of the polarization ensures that 
the NMR line of the nuclei probing the spins of the lat- 
tice is broadened on both sides, i.e. that the broadening 




FIG. 19: A gaussian staggered polarization of arbitrary inten- 
sity generated by a vacancy (center of the figure) on a 50 site 
kagome lattice. ^ is twice the spin-spin distance and the stag- 
gered vector is Q = (1/3, 1/3), i.e. the vector of the y/S x \/3 
kagome spin configuration. 
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is symmetric and scales with the apphed field, in agree- 
ment with our data. A non staggered polarization would 
generate only local fields pointing in the same direction 
and the NMR line would then be only broadened on one 
side, ie the broadening would be asymmetric (as in the 
orphan spin model). The actual shape of the broadening 
(gaussian, lorentzian etc.) depends on x'(f'), on the vec- 
tor Q and on the geometrical details of the nucleus-spin 
coupling. 

On the macroscopic stand, the vacancy generates a 
paramagnetic defect with a macroscopic susceptibility 

_ il-pUf,{def) 

where /J,gjj((ie/) is the effective moment of a defect. Al- 
though the exact nature of the defect is not well [estab- 
lished and is still the object of theoretical studies,E2l two 
scenarios may be conceived. 1) The defect is the sum of 
all the polarizations, in other words the magnetization 
M = whose susceptibility is paramagnetic. 

The 1 /T variation of the overall moment has then no in- 
tuitive explanation. 2) The vacancy generates a localized 
paramagnetic moment on the n.n. spins, which generates 
the staggered response x'(^) in the spin network. In 
this case, one needs to understand microscopically why a 
paramagnetic moment is created. 

Quantitatively, the expected NMR broadening is mag- 
netic and related to x* ^-nd ^ by 

^H^Ax*H^,l-p)^ (4) 

where A is the hyperfine coupling constant between the 
spin and the nucleus. ^ is a function of ^ and of the 
dilution (1 — p). The form of T depends on the form of 

The NMR width is sensitive to all variation in T of ^(T) 
and x*{T) in contrast with the macroscopic measure- 
ments. The surprising simple l/T variation of AiJ we 
observe by Ga-NMR in SCGO might indicate that once 
the correlations set in below 50 K, very little changes oc- 
cur for ^(T) and x*{T) due to the high frustration. Con- 
cerning ^(T), neutron measurements are very rewarding 
since they do infer that the spin-spin correlation is con- 
stant with temperature when T < 60 K (Sec. ^). 

At low dilution, where the interaction between the de- 
fects can be neglected, !F{1 ~ p) oc [1 — p). The evolu- 
tion of C^ef /NMR with the dilution reflects the one of 
JF(^, 1 — _p). Fig. [l^ indicates that Cdef/NMR is hnear 
with (1 — p) at low dilution, but deviates progressively 
from linearity at higher dilution. This deviation is indeed 
expected. It can be associated with an "interference" be- 
tween the polarizations induced by the nearby vacancies. 
This phenomenon, known now for several years, was first 
evidenced ift,the case of the RRKY polarization in the di- 
lute alloys.Eil It is very important to note that C^^eZ/jvAfi? 
likely extrapolates to when p ^ 1 or in the worst case to 
a value lower than the observed values of Cdef/NMR for 



p < 0.95. Since the SG state occurs at nearly the same 
temperature in all the samples and Tg varies the oppo- 
site way of 1 — p (bottom panel of Fig. p^ ), this clearly 
demonstrates that the origin of the SG state is not related 
to the dilution-induced defects. Whether a small residual 
p- independent Curie term associated with intrinsic (topo- 
logical) defects might explain the freezing is still a matter 
of speculation as higher Cr-concentrations than p > 0.95 
would be necessary to conclude. 

To summarize, our NMR and macroscopic suscepti- 
bility results agree perfectly with a picture of a defect 
built on the AF correlations of the frustrated network. 
In conclusion we now elaborate on the effective value of 
the paramagnetic moment detected in the low-T macro- 
scopic susceptibility. 

Following Eq.(||) and from the Cdef constant of Tab. 0, 
we deduce a nearly p-independent value for the defect ef- 
fective magnetic moment of fj,^j:j{def) « 1.5/i^ smaller 
than the expected 3.87 /i^ value for a free S = 3/2 
spin. The weak ratio Cdef/C of Tab. ^ is then sim- 
ply explained in this viewpoint since Cdef/C = (1 — 
p)^ilj:f{def )/pnlfj:{Cr^+) « 0.15(l-p)/p. We can either 
speculate that the geometric frustration might diminish 
the efficiency of the polarization induced by the vacancy, 
indication that the nature of the defect is complex or the 
paramagnetic defect can be modeled into a S* = 1/2 spin 
quantum state. Also, in the case where /Zgjj((ie/) results 
from the sum of the oscillating spin polarization there is 
no reason to find a 3/2 value. Further theoretical work 
in this direction is indeed required to model the defect 
and achieve a quantitative understanding of our data. 



VII. SUMMARY AND CONCLUDING 
REMARKS 

The SCGO kagome bi-layer compound was studied for 
the first time through a local probe over a wide variety 
of Cr-concentrations. New and original properties were 
unravelled for the archetype of geometrically frustrated 
compounds. 

Ga(4/)-NMR allowed us to observe the kagome bi- 
layer susceptibility X frusta otherwise non accessible by 
macroscopic measurements. A maximum in Xfmstr oc- 
curs at a temperature Tmax ~ O.IO, robust to a dilution 
as high as « 20%. No gaped feature is evidenced. The 
maximum signals the appearance of short ranged mag- 
netic correlations in the kagome bi-layer, as observed in 
all the AF d = 2 Heisenberg systems. Tmax is a new 
energy scale for SCGO, which sets a new constraint on 
the theoretical models. 

A close examination of the dilution effects by Ga(4/)- 
NMR allowed us to establish, on a firm experimental 
ground, that paramagnetic defects are present on the 
kagome bi-layer. The defects stem from the vacancies in 
the spin network, i.e. from the substitution of magnetic 
Cr^+ ions with non magnetic Ga^+ ions. The defects are 
responsible for the low-T Curie upturn of Xmacro- The 
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macroscopic susceptibility therefore does not probe an 
intrinsic property of the kagome bi-layer at low-T, but 
mainly a property related to the dilution. Interestingly, 
the SG transition is detected precisely in the low-T up- 
turn of Xmacro- there are experimental signatures 
suggest ing .that this transition is an intrinsic feature of 
SCGO,n3E3 this means that the defects do not trigger the 
freezing, but do freeze at T < Tg. The origin of the SG 
transition remains to be elucidated. 

Finally, we proposed a mechanism that justifies the ex- 
istence of the paramagnetic defects, which, to our knowl- 
edge, is the only one that gives a consistent interpretation 
for both the Ga(4/)-NMR and the macroscopic suscepti- 
bility results. Our NMR data point at a defect which pos- 
sesses strong analogies with the defects observed in most 
of the AF correlated systems — where the presence of a 
vacancy (or a magnetic impurity) generates a staggered 
response from the spin network. Up to now, the dilution 
in SCGO, or more generally in the frustrated systems, 
has been considered a parasitical effect to be minimized 
in order to measure the intrinsic properties related to the 
geometric frustration. In light of our results, it is interest- 
ing to reverse the problematic. The study of the response 
from the magnetic lattice to a non magnetic vacancy can 
be an indirect way of probing the intrinsic properties of 
the frustrated network such as £,{T) or x*{T)- The ex- 
perimental study and the theoretical models to describe 
this response are a new domain of investigation for the 
geometrically frustrated systems. 



where gs is the degeneracy of each level and j3 = l/ZcsT. 
The sum is taken over S" = 0, 1, 2, 3. ..2s. We now apply 
an external field H to calculate the susceptibility. The 
degeneracy is lifted by a Zeeman splitting E„ig (H) = 
g^gHms so that 



Z{H) 



S=0,2s ms = -S,S 



N 



where m^ = —2s, ...,0, ...,2s are the Zeeman levels for a 
spin S. From the expression of 2{H) we then derive the 
susceptibility of N isolated spin pairs 



E E mlexp{-(3Es) 



^9 H-B S=0,2srns = -S,S 



ksT E gsexp{-PEs) 

S=0,2s 



The susceptibility of a Cr(4f„i)-Cr(4/„i) spin pair (s 
3/2) per formula unit of SCGO is therefore 



Acknowledgments 

It is a pleasure to thank J. Bobroff, H. AUoul, A. Keren, 
C. Lhuillier and F. Mila for fruitful discussions. 



APPENDIX: SUSCEPTIBILITY OF AN 
ISOLATED SPIN PAIR 



X.pa 



x_ 

N 



s{s + l)kBT 

2exp{-(3J) + I0exp{-3PJ) + 28ea;p(-6/3J) 
1 + 3exp{-l3J) + bexp{-3(3J) + 7exp{-Qj3J) 



We consider pairs of spins s with an AF exchange 
constant J. Their states are labeled by their total spin 
quantum number, S = 0, 1,2, ...,2s, and lie at energies 
Es — J S{S +l)/2 relative to the ground state. The par- 
tition function for N independent pairs {2N being the 
total number of spins involved in pair formation) is 



Z = 



X! 9sexp{-l3Es) 



N 



where we have introduced the effective moment /^e// — 
\/9^t^%s{s + 1). The value of the AF exchange constant, 
i.e. of the singlet-triplet gap, is J = 18.6(1) meVJ23 To 
account for the dilution effects, Xpair weighted by p| j^. , 
the statistical probability of having a pair of Cr(4/^i)- 
Cr(4/i,i) per formula unit of SCGO. The value of Pif^- is 
known from neutron refinements (Fig. ^). Note that in 
the limit T 3> J/kB, we recover the sum of two s = 3/2 
paramagnetic susceptibilities: 2p\^ fj,^j:j:/3kBT. 
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